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1. INTRODUCTION 

Thn d o c -act describes a computer iiM«WM program which is used to estimate the effects of a proximate diffrac- 
ting fence — fee performance of parabufcnd anfeMas. 

The computer program is written in FORTRAN language for running on an IBM 3081 computer system at Goddard 
Space Plight Center, GreenbeH, MD. 

The physical problem, mathematical fonnufartk*. and coordinate reference* are described in the general description 
section. 

The main control structure of the p ro g ra m a ad fer on— of dr individual subroutines me discussed in the pro- 
gram description sectioo. 

The Job Control Language set-up and pragniat fetroett— ase p rovi ded in fee user’s i as tro c ti on section to help users 
to execute the present program. 

A sample problem with an appropriate output listing is made available as an illustration of fee usage of the program 

Finally, a summary and c o mprehe nsi ve results are picitatad in the last section. 

2. GENERAL DESCRIPTION 

The computer program is written to compute the effects of a nearby fence on the antenna receiving characteristics. 
Specifically, a purpose of this investigation is to find out how the mabdobe and the first sidelobe of the paraboloidal 
antenna’s reception pattern change due to the fields diffracted from the edge of a fence. This diffraction field could be 
described by n Somme rfeld's half-plane solution if the fence la not very far from the antenna and the angle from the 
antenna center to the bodrendi of the fence is large. (For details, see NASA TM-84996, R.F. Schmidt, “A Radio- 
Frequency Analysis of Paraboloidal Antenna Located Near Diffracting Fences.”) In order to estimate the effect of dif- 
fracting fields on the performance of an antenna, a focal plane reception pattern of the incoming field which includes 
an incident plane wave and a diffracted field from a fence is needed. This pattern shows how the mainlobe and 
sidelobes are distorted due Jo the introduction of the diffracting field. The pattern of a normally incident plane wave 
(without diffracting field) !* well-known by studying fee distrotion of the pattern of a Sommerfeld half-plane solution, 
which includes an incident plane wave and a diffracting field, the effect of the diffracted field from a fence on the per- 
formance of an antenna is estimated. 

The focal plane reception pattern are obtained as followed. An aperture field distribution is obtained by. calculation 
of a So— naffeld’s half-plane solution at the aperture region in the antenna reference frame. The Fourier transformation 
of fe» spart an fetid dferihatfea is fee inoal pfonr reception patt er n (Details, see appendix A of this documentation.) 

2.1 COORDINATE SYS’ITEM 

Figure 1 illustrates an inertial reference I tame (\.j./)|. antenna reference frame (x.v.z) A . ami fciue ivteien. • : • 
(x.y.rlp. The orientation of the antenna reference frame relative to the inertial reference frame is represented by rare . 
Eulcrian angles (oJJ.y) A and a translation vector T A . Similarly, the orientation of the fence reference frame relative to 
the inertial reference frame is represented by three Eulcrian angles (o,/?.y)f: and a translation vector Tp. 

Using these transformations, the desired quantities can be easily transformed among these three reference frames. 

The notation R ,- * f denotes the Eulcrian rotational transformation from inertial reference frame to fence reference frame. 
Similarly, the notation R ,- ’ A means the Eulerian rotational transformation from inertial reference frame to antenna 
reference frame, etc. Some useful relations among these transformations are cited below. 

The cascade transformation: R A “* f = R ,_ * F R A “* 1 . The inverse rotational transformation: R ,- * F = (R p “* *) transpose. 
The basic rotation transformation: 




cosy cosa - cos/? sina siny cosy sina + cos/? cosa siny 
-siny cose - cos/? sina cosy - siny sina + cos/? cosa cosy 


siny sin/? 
cosy sin/? 


sin/? sina 


- sin/? cosa 


cos /? 


Where (ojl.y) are :hm Fttlrrian rotational angles. 

The inertial rvl . •• • »*.»•>. n> (ho document is an earth-fixed coordinate system. 
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2.2 PLANE-WAVE 

The incident wave could he E-potarized, H-potnriaad or mirtw* of both polarization. 


E-polarization 

E * (-coao ua0, - sino sin/), co a/J) e _iK * (1) 

H «<- tine, coso, 0) e _iKS (2) 

H -polarization 

E = (sino, - coso, 0) e _iKS (3) 

H = (- coso sin/?, - sino tin/?, cos/J) e -iKS (4) 

The phase factor of the plane wave is 

e ~«KS m e ~3c • (* cut* cat# + y liw eo$fi * i unfit 


In our program, we use the incident angle (a, and 4 k pohrtuMi angtedi in the inertial reference frame. The 
incident plane wave is plane wave ** coadj (E-pel) + dsprntffiig on the polarization angle d(, the plane 

wave could be E-polarized, H-potamcd, or other lisnrtjr pdeM plane wave. 
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23 DIFFRACTION FIELD 

Sommerfeld’t solution fw » Am 
system unit as 

(E-plane pobriaAon) 

E, - -H*JM 
fiy - H,«n*i 


|4m*-w«w by a balf-pkalBa ls given in Ac Gaussian 


. (6) 

m 


-E 

^ , S-J-co.peattro^-t^fGW - QM 
V« 

-si ^ 

H . S_l e auK«*# * «*A. 

* /« 

{* . HW ♦ <H«01 f 'O* f } 

\v. -v ^;-y* 

■-si.' „ 

f U « «_t e SC (roo# - litaS). 

y V« ■: . 

{<”“ KK,) - «<l>] +i v'j 7^5 “ f “ | } 

H* * o 

The companion expressions for the H-plane polarization are S‘ vcn below as 
(H - plane polarization) 

-M 

c „ £__A » iK (r co*/» - z linf) . 

vr« 

{sin a [G(p) — G(q)l «*!”!} 


(*) 


(9) 


(10H 

••• — .v&ai 

(ID 


(12) 


c _ ~ c 4 iK (r co«0 - z *»n/j) . 

“y n/v 

{cos « I«p) + c«,)l +'x/£^ cos f f } 

r., - 0 

H* • »in(* 

H y *■• - E, stn? 

-B 

H , £_4. cos p e iK(rco ^ _ * ,i, ^ , lO(P) + 0(9)1 
1 \f* 


(i? 1 

< u 

(15) 

(16) 

.07) 
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Where G(«) - e~ *•* [/* J U (-«•>] > fgr« J JJj V** *. 

rotates the PonmJ mudml wfetead mmatmltt >1 

U*M - t ,«*»**«; 

* it , **m* -4 • 

«gMt) * 1'tdlNi't'M.' 

» 1 •UdmHMt# : 

•nd q - - ^$4?' 

F • - 

The rn d d wi atgtesf*. 0) ami dre cylindrical coomb nntw (r. f , r) wed in die above equations are in fence reference 
frame. 

Once do twgnrtc Acid comp mi 

(H x , H y , H,V * sea dr H«-h + ate dp Hn-poi (18) 


have been detoptead, Ike transformation of 

(H», Hy. HJfc * *** a (H,. Hy. HJf (19) 

is wed to obtain the magnetic field of the antenna aperture in die antenna reference fame, where dp is the polarization 
angle in the fenocre fcre no e frarne. • v '•#*,’ 4 ' 

2.4 TRANSFORMATION AMONG REFERENCE FRAMES 

In order to use Sommerfeld’s half-plane aolndoo, the plane-wave incident angles and polarization angle should 

transform from the inertial reference frame to the fence reference frame. 

* 

The incident angle pair ( o ,<0) corresponds to a unit vector as 

(x, y, z)i = (cosa cos/), sina cos/), sin/J)i (20) 

The angles of arrival of the plane-wave in the fence reference frame are obtained by the rotation transformation 

(x. y. z) F = (x. y. z) t (21) 

The translation transformation is ignored here since only angles arc of concern. The arrival angles ol the plane-wave in 
the fence reference' frame arc found by solving the equation below. 

(x. y, z)p * (cow cos/), sina cos/), sin/% (22) 

For the transformation of the polarization angle 4j from the inertial reference frame to the fence reference frame 
linearly-polarized magnetic field with polarization angle d| in the inertial reference frame has components 

(H x , Hy, H z )| = cosdi (- sina, cosa, 0)j + 

(23) 

sindi (— cosa sin/), - sina sin/3, cos/))i 


Those magnetic field components transform into fence reference frame by 

(H x . H v - R«~ F (H x , Hy. H*), (24) 

Thr poterioaim »»gle dy>. in fence re fer en c e frame, are found by wiving dw equation 
(H x , Hy, H z )f «> cnedf (- um, rasa, (% -f 

smdp ( - cma tmfi. - nina ■ ted, cotfih- ' (25) 
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2.5 FOCAL PLANE RECEPTION PATTERN - A ,.. 

The antenna analysis program calculates the focal plane reception pattern of the incoming diffracted electromagnetic 
field by Fourier transformation of a Somme rfeld half-plane solution (which is in antenna reference frame) at the antenna 
aperture region. The pattern given in decibels (dfi), is obtained by 

P - 20k>g| /Hp(e)H,da| (25) 

aperture 


Where H p (0) if « plane wave function and H, it fos jfMUfeft At id daftributio* which is approximated by a Soramerfeld 

hatfitoae aafctk*. 


3. PROGRAM DESCRIPTION 

This section describes a FORTRAN program which is used to calculate a focal plane reception pattern for a 
paraboloidal antenna and a nearby fence. 

The program hierarchy chart is shown in Figure I . This chart shows the flow of the program. The COMMON block 
cross-reference matrix is shown in Figure 2. This matrix shows foe COMMON blocks used in each subroutine. 

The subroutine functional descriptions and input parameters list are also included in this section. The program listing 
is provided in appendix B. 


3.1 SUBROUTINE FUNCTIONAL DESCRIPTION 

This section describes the functions performed by each subroutine. 4 

-- - 

MAIN The main routine controls the flow of the program. MAIN routine first calls subroutine ATIN to read the in- 
put parameters. MAIN routine then calls subroutine ATDEF to define the constant values in the program. 
The third subroutine called by MAIN is subroutine ATEXPL, which Is used to print out the significant 
parameters used id the program. The actual calculation is performed aftet these calls. 

ATIN This subroutine provides the input values for the program. 

ATDEF This subroutine defines constant values and converts physical units. 


ATEXPL This subroutine prints out sonic input parameters for the user s information and record. 

ATAPER This subroutine subdivides the reflector aperture into small differential areas and evaluates the coordinates, 
unit normal vector, and differential area for each small area at the aperture surface. 

ATROP This subroutine provides the rotation operators for the use in the rotational transformations among inertial, 
antenna, and fence reference frames. 


ANGROP This subroutine defines. Eulerian rotation operation A for the rotational transformation 


Where 


(cosy coscr — cos/? cosor siny) (cosy aina + cos fi com silty) (siny sin/?) 

( — s»my cma - co*0 sin o cosy) (- siny sina + cos /? cos a cosy) (cosy sin/?) 
(sin/J cosa) (— sin/? cosa) (cos/?) 
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INPUT: Eukrian angle* «, p, y. 

OUTPUT: Rotational Operator marttix A. 

TRANSP This subroutine transposes the rotation operator to obtain an inverse matrix for the inverse rotation 
transformation. 

This subroutine c o mpot es the inverse matrix 



CROSS This subroutine performs matrix multiplication. 
C (3,3) - A (3.3) B (3.3). 


Where 


C (ij) «* k< ^j 4 A (i,k) B (kj) ^ 

This subroutine performs rotational transformation £romi>ne reference frame to another reference frame. 




where A is Eulerian rotation operator. 

ATXPW This subroutine transforms plane-wave incident angles of arrival and polarization angle from inertial 
reference frame to fence reference frame. 

ATRAT This subroutine performs the necessary transformation to transform the cartesian coordinate point in antenna 
reference frame to the cylindrical coordinate point in die fence reference frame. 

This subroutine computes 

(x, y. z)j * [r a ~*] (x. y, z) A + T A (1 

(x. y, z) F * Ir , '* f ) (x, y, z)j - T F (2 


(r. 6 , *> F - (7 xp 2 + xp 2 , tan -1 


-l yp \ 


where Tp are translation vectors of the antenna reference frame and the fence reference frame to* the 
menial reference frame, respectively. The first equation transforms a cartesian coordinate point from the 
*r>tmna reference frame to the inertial reference frame. The second equation transforms that same point from 
rt*.* inertial reference frame to the fence reference frame. The third equation transforms that point from tht 
c^fttioar coordinate system to the cylindrical coordinate system (Sommerfeld solution equations arc 
devohtd h\ » <.v!indric*l coordinate system in the fence reference frame.) 
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ATPL 


ATCS 

ATCY 

ATRFA 


ATPWH 


ATCORR 


This subroutine computes G(p) and G(q) values in the Soramerfeld soiotion equations. 
This subroutine evaluates 

G(a) - e _ “ : [i/n e ~ U (-a) •* sfMa) J * & dp] 


wbwe UU) • |. % > fe 

3, « < fl 

*§•**) • i. a > 0 




and 


<| « - (2Kr **0)* cos^-~ 


: p * ■-;Wh Ms i 


a 2 c* 

0 VTsrt 


dt 


x e** 


The Fresnel integral f — — dt is evaluated at the subroutine ATCS, . 

, J ° V2nt 

This subroutine evaluates Fresnel integrals. 

C(X) * INTEGRAL (COS(T)/SQRT (2 • PI • T) SUMMED OVER T FROM 0 TO X) 
S(X) = INTEGRAL (SIN(T)/SQRT (2 • PI • T) SUMMED OVER T FROM 0 TO X) 


-4 

J 


This subroutine evaluates H -plane Sommerfeld solution equation*. 

This subroutine transforms the H -plane electromagnetic field components of the Sommerfeld solution from 
fence reference frame to antenna reference frame. 

This subroutine evaluates 


(H„ H,. H,)a - (R^ A ] (H„ Hy, H z ) f 

This subroutine evaluates the plane-wave functions H p at every point on the reflector aperture 

H p **•(— sina, coso, 0) e~ lKS (I) 

C “iKS = c ~*iK • (x cos* cosfl + y lino cos/? + z »in/J) (2) 

angles o, 0 describe the direction of the plane wave 

This subroutine caktibte* focal ph*c *vey*iot> p at t ern by the following equation 
Pattern * 20 log j j H P d$ J 

aperture 

where t% Sommerfeld solution and H p is plane wave hinmon 
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MAIN 


ATIN 


ATDEF ATEXPL ATROP 



ATXPW ATAPCR 

I 

ROT 


ATRAT 


ATPWM ATCORR 


ATCY 


ROT 


ATPL ATRFA 


ATCS ROT 


FJguri 2THI«rarchy OwrtT 
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ATIH ATBKDT ATDEF ATEXPL ATROP 
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3.2 INPUT PARAMETERS 

The names and description* of the program mput variable* are listed helow 


lnbl»- 2 Sa<tn*' <um* Descriptions of the Program Variable* 


Common 

Parameter 

Parameter 

! Description | 

4 -7* . .. . . 

Block 

Name 

Type 

INPUTl 

S1GMAM 


AnictvM iv Hector radius ] 


SIGMAO 

j 

Original value of sigma 


XIJ 

K'K 

Integration control constant 


XUM 

R*8 

Wavelength 


FR 

R*8 

Frequency in GHZ units 


F 

R*8 

Focal length of ideal paraboloid 

INPUT2 

ALPHA! 

R*8 

Plane-wave azimuthal angle in 
inertial system 


BETAI 

R*8 

Plane- wave polar angle^in iner- 
tial system "* J 

j 


DELTAI 

R*8 

a 

Plane-wave polarization angle in 
inertial system 


ALF1 

R*8 

Eulerian angle in antenna system 


BET1 

R*8 

Eulerian angle in antenna system 


GAM1 

R*8 

Eulerian angle in antenna system 

INPUT3 

ALF2 

R*8 ! 

; i 

Eulerian angle in fence system 


BET2 

R*K 

Eulerian angle in fence system 


GAM2 

R*8 

Eulerian angle in fence system 


XT A 

R*8 

Translation vector component for 
inertial-antenna correspondence 


YTA 

R*8 

i 

Translation vector component for 
inertial-antenna correspondence 


ZTA 

R*8 

Translation vector component for 
inertial-antenna correspondence 


XTF 

R*8 

Translation vector component for 
inertial -fence correspondence 


YTF 

R*8 

Translation vector component for 
inertial-fence correspondence 


7.TF 

R*8 

Translation vector component for 
inertial-fence correspondence j 

_ - — - - J 



INPUT 




Table * 


ANGA 

1*4 

angle o of plane wave H p 

ANGB 

R*8 

angle P of plane wave H p 

j 

DB 

R*8 

increment of angle p 

NANG 

1*4 

number of angles to be evaluated 
in the pattern 

XDF 

R*8 

If XDF * 1, the Sommerfeld 
solution equations arc fully 
calculated 


1 


If XDF the cylindrical 
diffraction part of Sommerfeld 
solution it discarded 
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4. USER INSTRUCTION 

This section comprise* the user's guide f«*» the Mtotvu analysis program. 


Presently, the program ctatt Mt tin fvna of rvxi executable loud modules in library XRHSC. ANTENNA. LOAD lor 
all hut the MAIN and BLOCK DAI A subprogram. Since these cannot he linked in automatically, only object modules 
wciv generated for them Thera a* located in libraiy XRHSC. LIB. OBJ, which also includes the object modules 
generated w the praams of c noting the members of the nonexecutable load library. 

Figure A shows the JCL to be used to recompile the program when no nonexecutable load library module is to be 
created. Figure B shows the JCL to be used when a Mcuthor r> to br created and added to the non-executaWc load 
library. Cataloged procedures in library 0*Yl clantMc die need for detailed copying of this JCL 

(members FORTVOOM and fORTVAWf* 

C 

For executing the pgUNi. W, * and ye y eu en dt tra . has been created. PAT is located in library 

XRHSC. LIB. ChiTL Ura* FAT ra etptrradnt to uttag die JCL in F«ura C. 


5. SAMPLE PROBLEM 

This section provides a sample problem to simulate the radiation pattern of a 9-meter paraboloidal reflector with F/D 
ratio of 0.44 with feed' at focal point. The incident electric field is **«t»med to be linearly polarized in the Z-dircction 
in the initial coordinate system. Operating frequency is 2.0 OHz with integration control constant XLI set to 0.5. The 
Eulerian angles in antenna coordinate system are (180.90,90) and in the fence coordinate system are (90,90,270). The 
translation vector compon e nts from inertial to aperture origin is (0,1000,0) and from inertial to fence origin is 
(0,2200,-700) in cm units. The radiation pattern is computed at aperture transmitted region with angular coordinate Q 
varying from 93.0 to 87.0 at 0.1 decrements. 

• 

All the default input values are set up in the BLOCK DATA of the program ANTENNA. FORT(PAT). The user 
may use the namelisted input in the data file ANTENNA . D ATA(CASEP0 1 ) to override the default input values. 

Computer print-out listings are shown in Figures D-G. The focal plane reception pattern is plotted in Figure H. 


6. SUMMARY 

A comprehensive test of the program has been performed. Figure 1 shows that the focal plane reception pattern with 
a fence oriented orthogonally to the paraboloid axis and located directly in front of the antenna )6 meters away with the 
antenna lower half blocked by the fence. Figures 2-10 show a series of focal plane reception patterns with the fence 
lower by 30 cm, 100 cm, 200 cm, 300 cm, 400 cm, 450 cm, 600 cm, 700 cm and 800 cm, respectively. It is clearly 
illustrated that the focal plane reception patterns are distorted for all the cases that the antenna arc blocked by the fence. 
The radius of the antenna is 450 cm. The distortion depends on how much are* is blocked. For example, Figure 6 
slight distortion of the reception pattern for a slightly blocked antenna, and Figure I shows great distortion of 
rcccjwn fur a bait blocked antenna h r ai*. shown thai all the unblocked antenna cases, the reception pattern 

are uikliMoned illustrated in figures 7 - 10. the first side lobe level is seen to be at —17.6 dB with respect to the 
main beam peak due to the fact that the effect ol space divergence between a focal point feed and the parabolic reflec- 
tor, and feed directivity, were not included in the program initially. 
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For the following two cases, although the antennas are not blocked by the fence geometrically, but the diffraction 
effects are shown on the sidelobe of receptkm pattern nevertheless. Figure 1 1 illustrates that with the fence 1 27 meters 
away from the antenna and 5 meter* lower from the level of center of antenna, with this orientation the angle of incom- 
ing diffraction field is 2.15 degrees and correspondingly the reception pattern has distortion at 6 •* 87.75 degrees. The 
same effect is also shown in the second case. In Figure 12, with the fence located at (430 m, -9 m), the angle of in- 
coming diffraction field from the fence is 1.15 degrees. The reception pattern cleariy shows distortion at 0 - 88.85 
degrees. 

In this investigation we found that as long as the antenna is not blocked by the fence, the main lobe of the reception 
pattern will not be distorted. However, the sidelobe will be distorted if the angle of the incoming diffraction field is 
roughly equal to the inclination angle of that sidelobe. If the angle of the incoming diffraction field is large compared 
to the inclination angle of the sidelobe concerned, the diffraction effect from the fence is not shown in the focal plane 
reception partem. 
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APPENDIX A 


RELATION OF FOCAL PLANE RECEPTION PATTERN AND 
FOURIER TRANSFORMATION Of THE FIELD DISTRIBUTION 
AT ANTENNA APOtnM REGION 

The aperture fidd liirt rfcp tiuu H, is approximated by a SommeffekJ half-plane aoiuttea described in the antenna 
. Reference frame Hu* agamva field distribution includes the incident plane wave and the incoming diffraction field from 
' the fence. In oa dt l an estimate the effect of the diffracting fence on the performance of the antenna, a focal plane 
reception pattern eosmponding to this aperture field distribution is needed to illustrate how the introduction of the dif- 
fraction field from the fence changes the focal plane reception pattern. 

In FignmA'l, P is an arbitrary point on the focal plane of antenna, A is the vertex of antenna, and F is the focal 
poin t . The angle between PA and FA is 8 . For a small angle approximation, a plane wave with an incident angle 8 
from llte opposite side will focus on point P after reflection from the paraboloidal antenna. Alternatively, the phase dif- 
ference Amction for any point on the aperture plane is e 1 *? * R for that corresponding point P on the focal plane. 

(Here kp is the unit wave vector of the plane wave with an incident angle 8 and R is a position vector to the-aper- 
ture plane.] The field at point P for the aperture field distribution H, can therefore be calculated by the integrate# of 
the aperture field distribution Hs multiplied with the phaae difference function e iK r ' R over the entire aperture region. 

/ H, • e*» • R da . ‘ / i \ ' 

■fCftUIC .^4 

Similarly, for any other point Q on the focal plane, there is a wave vector kg in the opposite side of AF with an 
incident angle equal to FAQ. 

The field at point Q is # 

/ H* - ‘ * da 

aperture 

Thr f.va! reception fwiurn fi*» ancmnv field divtribution H s is now represented as the Fourier transforma- 
tion o' H v ;ii the a|vrii:?f rvi'*- ; 

/ H s • e* ‘ R da 

aperture 



APPENDIX B 


The entire antenna analysis program is listed below. 
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KXKX TSO FOREGROUND HARDCOPY MKKX 

DSNAME *XRHSC . L I B . MVS . Cl I ST C FORTVCOM > 


OOOOOOIOPROC 2 SUFFIX MEMBER TERMON MIN<0) SEC<30) 

0000001 1IF 8SUFFIX « XXX THEN SET SUFFIX* 

000 00 020 ACCESS JC ACCT(SPONS) BOX(BOX) 

00000030QED TEMPJOB.CNTL NEW EMODE „ TTUC ,. MTU . 

00000040IN //4SYSUID.C0M JOB ( FH011, B22, 2) , ' UB4SUFFIX. (4MEMBER. ) * , TIME* (AMIN . ,4SEC. ) » 
00000050IN // MSGCLASS*U,CIASS*0,N0TIFY*4SYSUID 
00000060IN /KJOBPARM QUEUE«FETCH 

00000070IN /KPROCLI B*XRHSC . MVSPROC . CNTL , DDC , tv AUTKU ua 

OOOOOOSOIN // EXEC FORTVCOM, OUT* * X* , SUFFIX* f 45UFFIX. • ,MEMBER*8MEMBER,PREFIX*ANTENNA 
00000090IN // EXEC NT SO 
OOOOOIOOITERMON. SCHEDULE K 
00000110END N 


4 



xxxx TSO FOREGROUND HARDCOPY xxxx 

DSNAME=XRHSC. MVSPROC. CNTL C FORTVCOM) 


//FORTVCOM PROC USRID=XRHSC, PREFIX' ,OUT=’X» , BLKSIZE=7265, 

// TERMOUT= , x»,ERRLVL= , NOFIPS,FLAG(E)' 

//SOURCE EXEC PGM*FORTVS.REGION*2048K,COND*(4,LT), + 

// PARM* ' LC( 80 ) , 4ERRL VL , SOURCE , XREF , MAP , NODECK, NOL 1ST, OPT ( 3 ) * 

//STEPIIB DD DSN*SYS1 . FORTVS, DISP*SHR 

//SYSL1N DD DSN*4USRID. .4PREFIX4SUFFIX. . OBJC4MEMBER) , DISP=SHR, 

// UNIT=SYSDA, 

// DCB*( »RECFM*FB, LRECL*80, BLKSIZE s 3200) 

//SYSPRINT DD SYS0UT=40UT,DCB=(RECFM=VBA,LRECL*137,BLKSIZE*4BLKSIZE) 
//SYSPUNCH DD DUMMY, DCB*BLKSIZE-3440 

//SYSIN DD DSN*8USRID. .8PREFIX8SUFFIX. . F0RTC8MEMBER) , DISP'SMR 
//X SYSLIB DD DSN*4USRID. .8PREFIX&SUFFIX. .COMN.FORT ,DI$P*SHR 
//SYSTERM DD SY$0UT*4TERM0UT 


00000010 

00000020 

00000030 

00000040 

00000050 

00000080 

00000070 

00000080 

00000090 

00000100 

00000110 

00000120 

00000130 
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**** TSO FOREGROUND HARDCOPY M * * 

DSNAME«XRHSC . L I B . MVS . CL I ST < FORTVADD) 


OOOOOOIOPROC 2 SUFFIX MEMBER TERMON MIH(O) SEC(30> 

00000020IF ASUFFIX « XXX THEN SET SUFFIX* 

00000030SET CC » CC 

00000040IF ATERMON « TERMON THEN SET CC* L 
00000050ACCESSJC ACCT(SPONS) BOX(BOX) 

00000060QED TEMPJOB.CHTl NEW EMODE 

22222 070IN "‘SYSUID.CAD JOB CBSP0NS,B22,2), , LIB4SUFFIX.<»MEMBER.) , ,TIME*<«MIN. < 
22222 moin msoclass*u,ciass*o,notify*asysuid 

00000090IN /KJOBPARM QUEUE "FETCH 
000001 OOIN /NPROCIIB-XRHSC.MVSPROC.CHTL 

00000110IN // EXEC FORTVADD, NBLK"4 , OUT ■ ' ** * SUFFIX* 'ASUFFIX . ' »MEMBER*AM£MBER, 

00000120IN // PREFIX-ANTENNA 

00000130IN /✓ EXEC NTSO 

00000140V 

00000150BCC 

00000160ATERMON. SCHEDULE K 
00000170END N 


**** TSO FOREGROUND HARDCOPY **** 

DSNAME«XRHSC . MVSPROC . CNTL ' < FORTVADD ) 


//FORTVADD PROC USRID a XRHSC, PREFIX* , OUT * '*' , BLKSIZE=7265> 

// NBLK*40 , TERMOUT = ' X * , ERRL VL * * NOFIPS , FI. AG( E) ' 

//SOURCE EXEC PGM=FORTVS, REGI0N*2048K, COND* ( 4 , LT) , + 

// PARM* * LC( 80 ) , A ERRL VI , SOURCE, XREF, MAP, NODECK, NOLIST, OPT ( 3) ' 

//STEPLIB DD DSN*SYS1 . FORTVS, DISP-SHR 

//SYSLIN DD DSN-AUSRID. . APREFIXASUFFIX.- . OBJ(BMEMBER) , DISP-SHR,, 

// UNIT *SYSDA, 

// DCB*( , RECFM-FB, LRECL *80, BLKSIZE*3200) 

//SYSPRINT DD SYSOUT *AOUT , DCB=( RECFM-VBA, LRECL *137 , BLKSIZE*ABLKSIZE) 
//SYSPUNCH DD DUMMY, DCB*BLKSIZE*3440 

//SYSIN DD DSN=AUSRID. .APREFIXASUFFIX. . FORTC AMEMBER) , DISP*SHR 
//SYSTERM DD SYSOUT-ATERMOUT 

//LIBNAME EXEC PGM«LI BRYGN2, C0ND*(4, LT) , REGI0N-175K, PARM* • « • 

//SYSLIB DD DSNAME=SYS1 .DUMMY, DISP*SHR 

//SYSOUT DD DSNAME*A8LIBM0D, DISP*(NEW,PASSJ,UNIT*3350, 

// SPACE* < 3200 ,(ANBLK,80),,, ROUND) , DCB«(RECFM*FB, LRECL *80, BLKSIZE*3200 ) 
//SYSPRINT DD SYSOUT *8 OUT , DCB*(RECFM*FBA, LRECL *81 , BLKSIZE*7 290) 
//SYSPUNCH DD DUMMY , DCB*(RECFM*FB, LRECL* AO , BLKSIZE*7280 ) 

//SYSIN DD DSN*SUSRID. . 8PREFIXASUFFIX . . OBJ(AMEMBER) , DISP*SHR 
//LINK EXEC PGM*LINKEDIT,C0ND*(4, LT),REGION*150K, 

// PARM* ' LIST, MAP, NCAL ,SIZE=(132K,12K)' 

//SYSLMOD DD DSN=AUSRID. .APREFIXASUFFIX. . LOAD( AMEMBER) , DISP*SHR 
//SYSPRINT DD SYSOUT* AOUT , DCB*(RECFM*FBA, LRECL*121 , BLKSIZE*3509) 
//SYSUDUMP DD DUMMY 

//SYSUT1 DD UNIT*3350,SPACE«(CYL,C1,1)) 

//SYSLIN DD DSN*AALIBMOD, D1SP*( OLD, DELETE) 


00000010 

00000020 

00000030 

00000040 

00000050 

00000060 

00000070 

00000080 

00000090 

00000100 

00000110 

00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 

00000210 

00000220 

00000230 

00000240 

00000250 

00000260 

00000270 
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ASEC.), 



kkkx TSO FOREGROUND HARDCOPY NXXX 

DSNAME'XRHSC.LIB.CNTL (PAT ) 


//XRHSCP01 JOB (FH011.B22,5),PATDF,MSGCLASS«A,TIME«U,), 

// N0TIFY*XRHSC,CLASS*0 

//X CLAS5-0.N0 TAPE, A i DEFAULT, Ei EVENING, F« WEEKEND. 

XXXJOBPARM LINES"100, QUEUE"FETCH 
/xpROCLIB'XRHSC.MVSPROC.CNTL 

//STEP1 EXEC PAT,VOL«ANT01,FILE«1,SIZE»’20<48K,256KS 
// REGION. LINK«6000K, REGION. QO'SOOOK 
//LINK. OBJECT DD X 
INCLUDE OBJUB(PAT) 

INCLUDE NENLIN(ZETA8) _ 

//GO.FTOBFOOl DD DUMMY 

//XGO.FTOBFOOl DD DUMMY ... , 

//GO.DATA5 DD DSN'XRHSC. ANTENNA. DATA(CASEPOl) , DISP-SHR, LABEL»(, ,, IN) 
// EXEC NTSO,MODE*ALL 


OOOOOOIO 

00000020 

00000025 

00000050 

00000040 

00000050 

00000060 

00000070 

OOOOOOBO 

00000090 

00000100 

00000110 

00000120 

00000150 


XXXK TSO FOREGROUND HARDCOPY KXKX •• . 

DS NAME«XRHSC . MVSPROC . CNTL (PAT ) 


//ANTASIM PROC N8LK«40,OUT«'X*,TERMOUT«*X*,FILE«1, *LA 
// , SIZE=U28K,12K', OPTION* 

//LINK EXEC PGM«IEHL,REGIOH«150K,COND*(4,LT). 

// PARM**MAP, LIST , SIZE*(ASIZE. ) , 80PTI0N* 

//NEMLIN DD DSN»SYS2.NEWZETA,DISP»SHR 

//SYSLIB DD DSN'SYSl . VLNKMLIB, DISP*SHR XRDMS 01/06/BJ 

// DD DSN*$YS1 . VFORTLIB, DISP'SHR 

// DD DSN^XRHSC .ANTENNA .LOAD , DI SP*SHR 

// DD DSN*SYS2.WP1055,DISP«SHR 

// DD DSN*SYS2 . IMSLS, DISP'SHR 

// DD DSN*SYS1 . FORTSSP, DISP«SHR 

// DD DSN«SYS1.MVTFTLIB,DISP*SHR 

// DD DSN«SYS2. VFORTLIB, DISP*SHR XRDMS 01/06XSS 

XX DD DSN*SYS2 . FORTL IB, DISP*SHR 

xxOBJLIB DD DSN«XRMSC . ANTENNA .CBJ(PAT' .DISP s SHR, 

XX DCF* ( RECFMr Ff. > RJICL *SC.K KS:ZE = 32 00) 


00000010 

00000020 

00000030 

00000040 

00000050 

00000060 

00000065 

00000070 

00000080 

00000090 

00000095 

00000100 

OOOOOllO 

00000120 

00000130 

00000135 


xxSYSLMOD DD DSN*&81CD«GD(G3FC.> , UHlT'DISK, SPACE=(6144, (8NBLK,20,1)), 

// DISP*(,PASS) 

//SYSPRINT DD SYS0UT*80UT, DCB*(R£CFH*FBM, LRECL *121 , BLKSIZE'350,9) 
//SYSTERM DD SYS0UT*4TERM0UT 

xxSYSUTl DD UNIT* (DISK, SEP*< SYSLMOD, SYSPRINT) ) , SPACE*(CYL , (1 , 1 ) ) 
//SYSLIN DD DSN*S80BJM0D. DISP*(MGD, DELETE) , DCB=RECFM*FB, 

// UNIT*DISK,SPACE*CTRK, 0) 

// DD DDNAME=CBJECT 

//SYSUDUMP DD DUMMY 

//GO EXEC PGM=X. LINK. SYSLMOD, CDND=C4,LT),REGI0N*800K 

//FT05F001 DD DDNAME S DATA5 

//FT06F001 DD SYS0UT*40UT , DCB*(RECFM*VBA, LRECL*137 , BLKSIZE'l 61 , BUFNO'l ) 
//FT07F001 DD SYSOUT«B,DCB*(RECFM*FB,BLKSIZE*7280,LRECl«8D) 

//FT08F001 DD UNIT«(6250, , DEFER) , LABEL«( 8FILE, NL ), VOL «SER*IVOL , 

// DCB*(RECFM*FB, LRECL*64» BLKSIZE*1024, DEN*3) , 

// DISP*( NEW, KEEP) 


//FT10F001 DD SYS0UT*80UT , DCB*X. FT06F001 00000300 

//SYSUDUMP DD DUMMY .00000310 
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OF POOR QUALITY 
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kxx* TSO FOREGROUND HARDCOPY xkx* 

DSNAME«XRHSC. MEMO. TEXT (END ) 

The complete output file of our sample problem is shown as « 


NAMES AND DESCRIPTIONS OF THE PROGRAM VARIABLES 

SIGMAM « *50.000 XXX ANTENNA REFLECTOR RADIUS 

SIGMAO * 0.000 XXX INITIAL VALUE OF SIGMA 

XLI « 0.500 XXX INTEGRATION CONTROL CONSTANT 

XLM * 14.990 XXX WAVELENGTH 

FR * 2.000 XXX FRENQUENCY IN GHZ UNJTS 

F = 500.0 XXX FOCAL \ ENGTH OF IDEAL PARABOLOID 

PLANE-WAVE ANGLES OF ARRIVAL IN INITIAL SYSTEM 

ALPHAI = J90.000 XXX AZIMUTHAL ANGLE 

BETA I = O.pOO XXX RIGHT ANGLE - POLAR ANGLE 

POLARIZATION ANGLE 

DELTAI « 90.00C 

EULEFIAN AN'OlEf PC® KOTATJON 

IN ANTENNA L I" 

AL FI = 180.000 BET 1 = 90.000 GAM1 = 90.000 

IN FENCE SYSTEM 

ALF2 = 90.000 BET2 « 90.000 GAM2 = 270.000 

TRANSLATION VECTOR COMPONENTS 

FROM INITIAL TO ANTENNA SYSTEM 

XT A * 0.0 YTA =1000.0 ZTA = 0.0 

FROM INITIAL TO FENCE SYSTEM 

XTF= 0.000000E+00 YTF = 0.220000E+04 ZTF « -.7000COE*05 


20 



vfct&iAL PAGE fS 
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INDEX 

e 

11992 





PB « 

93.000 

VO C « 

0 . 183998 E*05 

VOCL 

s 

0 . 852703E+02 

PB * 

92.900 

VOC « 

0 . 38 1938E+05 

VOCl 

* 

0 916398E+02 

PB « 

92. BOO 

VOC « 

0 . 558975E+05 

VOCL 

s 

0 . 999979E+02 

PB « 

92.700 

VOC * 

0.692831E+05 

VOCL 

X 

0 . 968 125E-*C2 

PB * 

92.600 

VOC « 

0.763899E+05 

VOCL 

X 

0.976607E+02 

-PB « 

92.500 

VOC « 

0.756856E+05 

VOCL 

9 

0.975803E+02 

PB * 

92.900 

VOC » 

0.6627C3E+03 

VOCl 

9 

0 .969269E+02 

PB *• 

92.300 

VOC * 

0 . 980307E+05 

VOCL 

9 

0 . 936309E+02 

PB • 

92.200 

VOC « 

0.217958E+05 

VOCL 

m 

0.867975E+02 

PB » 

92.100 

VOC « 

0 . 1 08810E+05 

VOCL 

X 

0.807339E+02 

PB « 

92.000 

VOC « 

0.972636E+05 

VOCL 

X 

0 . 939905E+02 

PB « 

91.900 

VOC * 

0.890626E+05 

VOCL 

X 

0.989921E+02 

PB = 

91.800 

VOC « 

0 . 117380E+06 

VOCL 

9 

0.101392E+03 

PB * 

91.700' 

VOC « 

0 . 193021E+06 

VOCL 

X 

0.103108E+03 

PB « 

91.600 

VOC * 

0 . 156805E+06 

VOCL 

9 

0.103907E+03 

PB * 

91.500 

VOC « 

0 . 159890E+06 

VOCL 

X 

0 . 103S00E+03 

PB * 

91.900 

VOC « 

0.139109E+06 

VOCL 

X 

C.102599E+C3 

PB - 

91.300 

VOC = 

0.922396E+05 

VOCL 

X 

0 . 992989E+02 

PB «= 

91.200 

VOC = 

0.283026E+05 

VOCL 

= 

0.890365E+02 

PB = ■ 

91.100 

VOC = 

0 . 57 3251E+05 

VOCL 

X 

0.951669E+02 

PB = 

. 91.000 

VOC = 

0 . ]6280?E*0£ 

VOCL 

= 

0.109233E+03 

PB = 

9 0 9 C 0 

VOC - 

0 2BC»e2 E+06 

VOCL 

= 

0 . 1 09093E+03 

FB - 

9 0 , £ 0 C 

VOC = 

Q . 9 18997 E+06 

VOCL 

= 

0 .112999E+03 

PB = 

90.700 

VOC = 

0.559590E+06 

VOCL 

s 

0 .1199576+03 

PB = 

90.600 

VOC = 

0 . 7 00155E+06 

VOCL 

c 

0.116909E+03 

PB = 

90.500 

VOC s 

0.839125E+06 

VOCL 

= 

0 . 118925E+03 

PB = 

90.400 

VOC = 

0.959788E+06 

VOCL 

s 

0 . 119598E+C3 

PB = 

90.300 

VOC <= 

0 . 105597E+07 

VOCL 

= 

0 . 120973E+03 

PB * 

90.200 

VOC * 

0.113237E+07 

VOCL 

s 

0 . 121 080E+03 
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PB 

S 

90.100 

VOC 

* 0.117996E+07 

VOCl 

s 

0.121657E+05 

PB 

c 

90.000 

VOC 

* 0 .119615E+07 

VOCL 

* 

0.121556E+03 

PB 

s 

89.900 

VOC 

« 0 . 118006E+07 

VOCl 

x 

0 . 121658E+05 

PB 

c 

89.800 

VOC 

« 0 . 115259E+07 

VOCL 

* 

0. 121081E+C3 

PB 

ft 

89.700 

VOC 

« 0.105627E+07 

VOCl 

s 

0. 120676E+05 

PB 

* 

89.600 

VOC 

« 0 . 955166E+06 

VOCl 

s 

0.11 96 01 E+05 

PB 

s 

89.500 

VOC 

> 0.856698E+06 

VOCl 

r 

0.118629E+05 

PB 

c 

89.600 

VOC 

« 0 .700509C+06 

VOCl 

* 

0 . 1 16908E+05 

PB 

• 

X 

89.500 

VOC 

• 0 . 5598*1 [+06 

VOCl 

ft 

0. 116961E405 

PB 

X 

89.200 

VOC 

* 0 <.19217 £+06 

VOCl 

s 

0.112669E+05 

PB 

S 

89 ICO 

VOC 

» 0 . 2850008*06 

VOCl 

x 

0 . 109097E+05 

PB 

c 

89.0CO 

VOC 

* 0 . 162816E+06 

VOCl 

c 

0.106256E+05 

PB 

c 

88.900 

VOC 

' 0.572262E+05 

VOCl 

X 

0 . 951516E+02 

PB 

x 

88.800 

VOC 

* 0.286991E+05 

VOCl 

X 

0.890966E+02 

PB 

c 

88.700 

VOC 

* 0.925082E+05 

VOCl 

s 

0 . 995256E+02 

PB 

= 

88.600 

VOC 

« 0 . 156616E+06 

VOCl 

x 

0 . 1 02569E+05 

PB 

c 

88.500 

VOC 

* 0 .155207 E+06 

VOCl 

r 

0.105818E+03 

PB 

r 

88.600 

VOC 

* 0 . 157096E+06 

VOCl 

= 

0 . 1 03923E+03 

PB 

s 

• 

88.500 

VOC 

« 0.165250E+06 

VOCl 

= 

0.103122E+03 

PB 


.88.200 

VOC 

= 0.117525E+06 

VOCl 

= 

0 .1 01602E+03 

PB 


88.100 

VOC 

= 0.861065E+05 

VOCl 

r 

0 . 986966E+02 

PB 

= 

88.000 

VOC 

* 0 . 671980E+05 

VOCl 

X 

0 . 956785E+02 

PB 

= 

87.900 

VOC 

= 0. 207148E+C5 

VOCl 

= 

0 . 805997 E+02 

PB 

= 

87 .800 

VOC 

* 0.2 19976 E+ 05 

VOCL 

= 

0 . 868675E+02 

PB 

= 

87.700 

VOC 

= 0 . 683606E+05 

VOCl 

r 

0 . 936862E+02 

PB 

c 

87.600 

VOC 

* 0.666055E+05 

VOCl 

= 

0 . 9667 02E+02 

PB 

r 

87.500 

VOC 

* 0.760108E+05 

VOCl 

s 

0.97617 5E+02 

PB 

- 

87.600 

VOC 

* 0.766696E+05 

VOCl 

X 

0.976925E+02 

PB 

s 

87.500 

VOC 

= 0 . 696880E+05 

VOCl 

X 

0.968382E+02 

PB 

t 

87.200 

VOC 

* 0 . 56 0065E+05 

VOCl 

c 

0 . 9696 66E+02 
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87.100 VOC 
87.000 VOC 


0 . 381897E+05 VOCL « 

0 . 182295E+05 VOCL * 


0 . 9163S9E+02 
0 . 8521 54E+02 
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T f t I "P'T T 



: W$- ;fWf « « 


YTF = 2200.0 
ZTF = -600.0 



0 DEGREES 





oooooooooo 
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***** TSO FOREGROUND HARDCOPY XXX X 
DSNAME*XRHSC. ANTENNA. FORT 

CXXXXXXXKXXXXXXKXXXKMXXKKKKKKKKKKKKXXKXKKKXKI 


(ATMAIN ) 


(V XKKKKNKKKKKXKKXKKKXXKKKOO 00001 0 

00000020 
00000050 
00000040 
00000050 
00000060 
00000070 

ooooooso 

00000090 

MAIN PROGRAM. • TO CONTROL THE LOOICAL FLOW OF THE PROGRAM 00000100 

00000110 

CXXXXXXXXXXXXXXXXXXXKKMXHXKKKVXXXKXXHXXHKXKKMXXMXXXKKXKKKKXXXKXKXXXXXKKK00000120 


PROGRAM OF ANTENNA-FENCE SIMULATION 
TASK. NO. 511 
MARCH 27, 1986 

GSFC ATR - RICHARD F. SCHMIDT 
SAR TASK LEADER - DR. MACHAEL KAO 
SAR TASK PERSONNEL - HWAI-SOON CHENG 


IMPLICIT REALX8 <A-H,0-Z> 

REALX8 KIX.KIY,KIZ,*FX,KFY,KFZ,KMX,KHY,KNZ 

COMMON /INPUT1/ SIONAM,SIOMAO,XLI,FR,F 

COMMON /INPUTS/ ALPHAI,BETAI,DELTA1,ALF1,BET1,GAM1 

COMMON /INPUTS/ ALF2, BET2,GAM2,XTA, YTA,ZTA,XTF,YTF,ZTF 

COMMON /INPUTX/ ANGA , ANGB , DB , NANG , XDF 

COMMON /ANTO/ PIE,SPM,RAD,XK,XIM,PA,PB,PDB 

COMMON /ANTROP/ RIA(S,37,RIF(3,S),RFA(3,3),RAI(S,3),RFI(3,S) 

COMMON /ANTXP1/ KIX, KIY, KIZ, KFX, KFY, KFZ, KHX, KHY , KHZ 

COMMON /ANTXP2/ HIX,HIY,HIZ 

COMMON /ANTXP3/ ALPHA, BETA, SINALP,COSALP,SINBET,COSBET 
COMMON /ANTXP4/ SINDEL ,COSDEL 

COMMON /ANTPA1/ XL C 11993), YL (11993), ZL( 11993), INDEX 

COMMON /ANTPA2/ XNI(11993),XNJ(11993),XNK(11993),DSS(11993) 

COMMON /ANTRAT/ RDPR,ZDPR, THETA 

COMMON /ANTPL/ GPR,OPI,OOR,ORI 

COMMON /ANTCY/ FKXR, FKXI »FHYR» FKYI « FHZR, FHZI 

COMMON /ANTRFA/ AHXR, AHXI, AHYR, AHYI , AHZR, AHZI 

COMMON /ANTRF1/ HXRU199S),HYR(11993),HZR(11993) 

COMMON /ANTRF2/ HXIC1199S),HYIC11993),HZI(11993) 

COMMON /ANTPH1/ HXAR(11995),HXAI(11993),HYAR(11993) 

COMMON /ANTPW2/ HYAIC11993),HZAR(11993),HZAI<11993). 

COMMON /ANTLP1/ L00P1 


CALL 

CALL 

CALL 

CALL 

CALL 

CALL 


ATIN 
ATDEF 
ATEXPL 
ATROP 
ATXPW 
ATAPER 


DO 1000 LOOP1 * 1, INDEX 
CALL ATRAT 
CALL ATPL 
CALL ATCY 
CALL ATRFA 


1000 CONTINUE 

DO 2000 L00P2 « 
CALL ATPWH 
CALL ATCORR 
PB * PB - PDB 
2000 CONTINUE 
STOP 
END 


1 , NANG 


00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 

00000210 

00000220 

00000230 

00000240 

00000250 

00000260 

00000270 

00000280 

00000290 

00000300 

00000310 

00000320 

00000330 

00000340 

00000350 

00000360. 

00000370 

00000380 

00000390 

00000400 

00000410 

00000430 

00000440 

00000450 

00000460 

00000470 

00000480 

00000490 

00000500 

00000510 

00000520 

00000530 

00000540 

00000550 
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non 


***** TSO FOREGROUND HARDCOPY **** 
DSNAME«XRHSC .ANTENNA . FORT 


(ATBKDT ) 


BLOCK DATA 00000010 

CXKXKXXKKXXKKKKNKXXXXXXXXXKXXKXKXKXKKKKXXXKXXXXXXXXXXXXXXXXXKKKXKKXKXXXX00000020 


BLOCK DATA « TO SET UP ALL THE DEFAULT VALUES 

CXXXKXXXXXKXKKKXXKKNXXKXXXXKXKKXXXXXXXXKXXMXXKXXKXXi 


00000030 

00000040 

00000050 

KKXXXXKXKXXKXXXXO 0000060 


IMPLICIT REAL** <A-M,0-Z> 

COMMON /INPUT1/ $IGMAM,3IGMA0.Xll,Fft,F 
COMMON /INPUT2/ ALPHAI , BETAI , DELTA1 , ALF1 , BET1 ,0AM] 
COMMON /INPUT 3/ ALFZ, BET2,0AM2,XTA, YTA,ZTA,XTF, YTF.ZTF 
COMMON /INPUTX/ ANGA,ANGB,DB,NANG,XDF 
DATA SIGMAM/450 . 0/ , SIGMA 0/0 , 0/ 

DATA XLI/0.5/,FR/2.0/,F/500.0/ 

DATA ALPHAI/90 . 0/ , BETAI/O . 0/ , DELTAI/90 . 0/ 

DATA ALF1/1B0 . 0/. BET1/90 . 0/.GAM1/90 . 0/ 

DATA ALF2/90 . 0 / , BETZ/90 . 0/,GAM2/27 0 . 0/ 

DATA XT A/ 0 . 0/ , YT A/ 1 0 0 0 . 0/ , ZTA/ 0 . 0/ 

DAT A XT F/ 0 . 0/ , YT F/ 22 0 0 . 0/ , 2T F/-7 0 0 . 0/ 

DATA ANGA/O . 0/ , ANGB/93 . 0/, DB/O . 1/ , NANG/61/, XDF/1 . 0/ 

END 


00000070 

00000080 

00000090 

00000100 

00000110 

00000120 

00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 



nnnnn 


xxxx TSO FOREGROUND HARDCOPY XXXX 
DSNAME-XRHSC . ANTENNA . FORT 


(ATIN 


) 


SUBROUTINE ATIN 00000010 

CKXKKKKHKKXKKKKKKKKKKXKKKKKXXKKXKHKKKKKKXKHKKXKXXKKKXKKKNKXXKKNHHKKXMKXKOO 000020 

00000030 

SUBROUTINE ATIN • PROVIDING THE INPUT INFORMATION TO SPECIFY 00000040 

THE GEOMETRICAL AND PHYSICAL REQUIREMENTS 00000050 

FOR THE PROGRAM 000000*0 

00000070 

CXXKKKXXXXXXXKKXKXXXKKXXXXKKKXKXKKKKXKKXXKXKKXKKKXKKKXXKKKXKXXXXKKXKKKKXOOOOOOBO 
IMPLICIT REALMS (A-H,0-Z) 00000090 

COMMON /INPUT1/ SIGMAM,SIGMAO.XLI,FK,F 00000100 

COMMON /INPUT2/ ALPHAI, BETAI , DELTAI , ALF1 , BET1 ,QAK1 00000110 

COMMON /INPUTS/ ALF2* BET2,GAM2 * XTA, YTA.ZTA, XTF, YTF , ZTF 00000120 

COMMON /INPUTX/ ANGA, ANGB, DB, NANG,XDF 000001S0 

NAMELIST /INPUT/ SIGMAM, SIGMAO, XLI , FR, F, 00000140 

1 ALPHAI, BETAI, DELTAI, 00000150 

2 ALF1, BET1,GAM1, ALF2, BET2,GAM2, 00000160 

3 XTA,YTA,ZTA,XTF,YTF,ZTF, 00000170 

4 ANGA , ANGB , DB , NANG , XDF 00000180 

WRITE (6,70) 00000190 

READ C 5, INPUT, END*30, ERR*40) 00000200 

30 WRITE (6,80) 00000210 

GO TO *0 00000220 

40 WRITE (6,90) 00000230 

*0 RETURN 00000240 

70 FORMAT C1H///1H^»READ IN PROGRAM INITIAL VALUES') 00000250 

80 FORMAT ('XXX END OF PROGRAM INPUT DATA XXX') 00000260 

90 FORMAT ('XXX WARNING ... INCORRECT INPUT !!! •) 00000270 

END 00000280 


4 


34 


oooo 


XKXX TSO FOREGROUND HARDCOPY mkmk 
DSNAME*XRHSC. ANTENNA . FORT 


CATDEF ) 


SUBROUTINE ATDEF 
CXKXXXXXXKXXKXKXKMKXKMI 


IKXKXXKKKI 


SUBROUTINE ATDEF . TO DEFINE CONSTANT VALUE AND CONVERT PHYSICAL 

UNITS 


IMPLICIT REALX8 (A-H,0-Z> 

COMMON /INPUT1/ SIGMAM, SIGMAO,XLI , FR, F 

COMMON /INPUT2/ ALPHAI , BETAI . DELTAI . ALF1 , BET1 ,GAM1 

COMMON /INPUTS/ AL F2, BET2 . GAM2 . XTA, YTA . ZTA , XTF , YTF.ZTF 

COMMON /INPUTX/ ANGA, ANGB . DB. NANO. XDF 

COMMON /ANTO/ PI E, SPM, RAD, XK, XLM. PA , PB , PDB 

SPM * 2 . 997925D10 

XLM « SPM/(FRX1.0D9) 

PIE ■ 3 . 1415926556 
RAD * PIE/180.0 
XK « C2.0XPIE)/XLM 
ALPHAI * ALPHAIXRAD 
BETAI = BETAIXRAD 
DELTAI * DELTAIXRAD 
ALF1 «= ALF1XRAD 


00000010 
KKX00000020 
00000030 
l 00000040 
00000050 
0000006 0 
■■■00000070 
00600080 
00000090 
00000100 
00000110 
00000120 
00000130 
00000140 
00000150 
00000160 
00000170 
00000180 
00000190 
00000200 
00000210 
00000220 
00000230 


BET1 * BET1XRAD 
GAM1 * GAM1XRAD 
ALF2 « ALF2XRAD 
BETZ * BET2XRAD 
GAM2 « GAM2XRAD 
PA « ANGAXRAD 
PB * ANGBXRAD 
PDB * DBXRAD 
RETURN 


00000240 

00000250 

00000260 

00000270 

00000280 

00000290 

00000300 

00000310 

00000520 


END 


00000330 


TS 



ORIGINAL PAGE IS 
OF POOR QUALITY 


**** TSO FOREGROUND HARDCOPY xxxx 

DSNAME*XRHSC. ANTENNA. FORT (ATEXPl ) 


SUBROUTINE ATEXPL 00000010 

CKXXXXXXXXXXXXXXXKXXXXXKXXXXXKXXXXXXXXXXKXXXXXXX»XXXKKXXXXKXXXX*«*X«KKKK00000020 
c SUBROUTINE ATEXPL . TO PRINT OUT SIGNIFICANT INPUT PARAMETERS 00000030 
C FOR THE USER'S INFORMATION AND RECORDS 00000040 

C 00000050 

CXXXXXXXXXXXXXXXXXXXXXXXXKXKXXXXKXKkKXK«X*KXX*XXftKXXX*XKXXXKXXXXXXXKXXXK00000060 
IMPLICIT REALX8 (A-H.O-Z) 00000070 

COMMON /INPUT1 / SIGMAM, S10KAO, XL I , FR , f 00000080 

COMMON /INPUT 2/ ALPMAI . Bf TAJ . DElTAI , At F] . K7 1 , GAMi 00000090 

COMMON /INPUTS/ AL F2 . BET?, GAM2, XTA, YTA, Z7 A . XTf . YT F , 2! F 00000100 

COMMON /ANT 0/ PI E, SPK, RAD, XK, XLM, PA, PB , PDB 00000110 

common / input x/ anga,angb,db,nang,xdf 00000120 

WRITE (6,10) CC000130 

10 FORMAT ( ‘2 M5X, 'NAMES AND DESCRIPTIONS OF THE PROGRAM VARIABLES ') 00000 140 
WRITE (6,20) SIGMAM 00000150 

20 FORMAT ( ’0',15X, 'SIGMAM « ',F8.3,' XXX ANTENNA REFLECTOR RADIUS ' 00000160 
1) 00000170 

WRITE (6,30) SIGMAO 00000180 

30 FORMAT ( 'O', 15X, 'SIGMAO « ',F8.3,' XXX INITIAL VALUE OF SIGMA') 00000190 

00000200 

XXX INTEGRATION CONTROL C0NSTANT00000210 

00000220 
00000230 

XXX WAVELENGTH') 00000240 

00000250 

XXX FRENQUENCY IN GHZ UNITS') 00000260 

00000270 

XXX FOCAL LENGTH OF IDEAL PARABOL000000280 

00000290 
00000300 


WRITE (6,40) XLI 

40 FORMAT ( ' 0 • , 15X, 'XLI « ',F8.3,' 

1 ') 

WRITE (6,50) XLM 

50 FORMAT ( • 0 ' , 15X, 'XLM « ',F*.S,' 

WRITE (6,60)FR 

60 FORMAT ('0',15X,'FR * *,F8.3,' 

WRITE (6,70) F 

70 FORMAT C0',15X,'F « ',F6.1,' 

1ID' ) 

WRITE (6,80) 

80 FORMAT (• 0 ',15X, 'PLANE-WAVE ANGLES OF ARRIVAL IN INITIAL SYSTEM') 00000310 


• , F8 .3, ' XXX AZIMUTHAL ANGLE') 

XXX RIGHT ANGLE - P-GLAR ANGLE') 


WRITE (6,90) ALPHAI/RAD 
90 FORMAT ('O' , 15X, ' ALPHAI 
WRITE (6,100)* BETAI/RAD 
100 FORMAT ( 15X, ' BETA I « ',F8.3,' 

WRITE (6,105) 

105 FORMAT (' 0 ', 15X, 'POLARIZATION ANGLE') 

WRITE (6,110) DELTAI/RAD 
110 FORMAT ('O', 15X, 'DELTAI * ',F8.3) 

WRITE (6,120) 

120 FORMAT ( 'O' , 1 5X , ' EUl ERI AN ANGLES FOR ROTATION') 
WRITE (6,130) 


00000320 

00000330 

00000340 

00000350 

00000360 

00000370 

00000380 

00000390 

00000400 

00000410 

00000420 


130 FORMAT < ' 0 • , 1 5 X , 'IN ANTENNA SYSTEM') 

WRITE (6,140) ALF 1/RAD, BET1/RAD, GAMI/RAD 

140 FORMAT ( ' 0 » , 15X, ' AL FI = ',F8.3,' BET1 = ',F8.3,' GAMI =',F£.3) 
WRITE (6,150) 

150 FORMAT ('0',15X,'IN FENCE SYSTEM') 

WRITE (6,160 ALF2/RAD.BET2/RAD,GAM2/RAD 

160 FORMAT ( • O' , 15X, 'ALF2 = ’,F8.3,’ BET2 = *,F8.3,’ GAM2 = ',Ff.3) 
WRITE (6,170 

170 FORMAT ( *0* , 15X, 'TRANSLATION VECTOR COMPONENTS') 

WRITE ( 6 , 1 £ O 

UC FORMAT ( ' 0 ' , 1 5X , ' FROM INITIAL TO ANTENNA SYSTEM') 

WRITE (6,190 XT A . YT A . 7.T A 

190 FORMAT ('0M5X,-XTA - 1 , F6 . 1 . • YTA * • , F6 . 1 , ' ZTA = 

WRITE (6,200) 


C0000430 

00000440 

00000450 

00000460 

00000470 

CC0C0480 

0CCCC490 

cooocsco 

00000510 
00000520 
C0000530 
00000540 
00000550 
0000056 C 
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OF POOR QUALITY 


200 FORMAT ( » 0 ' , 15X, • FROM INITIAL TO FENCE SYSTEM') 00000570 

WRITE (6,210) XTF,YTF,ZTF 00000500 

210 FORMAT ( ' 0 • , 15X, 'XTF« *,E12.6,' YTF » \E12.6,' ZTF * *,E12.6> 00000590 

WRITE (6,250) ANGA, ANGB,DB 00000600 

230 FORMAT ( • 0 • , 15X, ' ANGA • ',FB.3.' ANGB « ',F8.3,' DB = ',FB.3) 00000610 

WRITE (6 , 240 ) NANG, XDF 00000620 

2A0 FORMAT ( ' 0 • , 15X, ’ NANG * »,I5,' XDF « '.F8.A) 00000630 

RETURN 00000660 

END 00000650 


i 


v? 



uuou 


XXXX TSO FOREGROUND HARDCOPY XXXX 
DSNAME«XRHSC. ANTENNA . FORT 


(ATROP ) 


SUBROUTINE ATROP 00000010 

CXXXKXKXXKXKKKKXKKKKKKKKKKXMKXXKKKKKKKXKXKKKKKKKKKK*KKKKXKXKXXXXXXXXXKXK00000020 
SUBROUTINE ATROP . TO PROVIDE THE ROTATION OPERATOR FOR THE 00000050 

USE IN ANTENNA-INITIAL, I NI TI AL -FENCE OOOOOOAO 

SYSTEM, AND VICE VERSA 00000050 

„ 00000060 

CXKXXXXKKKXXXXKXXXXX'«K».«KXXXKXXXXXXX* XXX MX * ****** ♦»«***********«*X*XXX*X00000070 
IMPLICIT REAL *8 (A-H.O-Z) 00000080 

COMMON /INPUT 2/ AL PHAI , BETA1 , DELTAI , ALF1 , BET 1 , CAM1 00000090 

COMMON /INPUTS/ ALF.2. BET2, GAM2,XTA , YTA , ZI A, XTF . YTF . ZTF OOOOOIOO 

COMMON /ANTRCP/ RI A ( 3 . 3 ) , RI F ( 3 , 3) , RFA( 3, 5 1 , RAI ( 3, S ) . RFI ( 3 , S ) 00000110 

CALL ANGROP ( Al FI , BET! , GAM1 , RAI ) 00000120 

CALL ANGROP (ALF2, BE72.GAM2. RIF) 00000130 

CALL TRANSP (RAI , RIA) 000001A0 

CALL TRANSP (RIF, RFI) 00000150 

CALL CROSS <RIA,RFI,RFA) 00000160 

RETURN 00000170 

END 00000180 


OfflGfflW- PAGE » 
OP POOR QUALITY 


3X 


noon 


XKXX TSO FOREGROUND HARDCOPY xxxx 
DSNAME=XRHSC. ANTENNA . FORT 


(ANGROP > 


SUBROUTINE ANGROPLA, B,C, T1 ) 00000010 

CXXXXXXXXXXXKKXXM«XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKKKXXXXXXKX»»X«N>r « » * »x* o 0 0 0 0 020 

00000030 

SUBROUTINE ANGROP i TO DEFINE THE ROTATION OPERATOR FOR 000000A0 

THE GEOMETRICAL TRAN' FORMAT J ON 00000050 

00000060 

CXKKKKXKKXXXXXXKXXXXKXXKXXXXXXMX XX» («■•••*•*■••• *«»« XXXXXXXXXXXO 0 0 0 0 07 0 
IMPLICIT REALMS (A-H,0-Z> 00000080 

DIMENSION TK3.3) 00000090 

SA * DSIN(A) 00000100 

CA « DCOS(A) 00000110 

SB * DSI N( B > 00000120 

CB *■ DCOS(B) 00000130 

SC « DSIN(C) 00000140 

CC » DCOS(C) 00000150 

Tl(l,l) * CCXCA-CBXSAMSC 00000160 

Tl(l,2) « CCXSA+CBXCAXSC 00000170 

T1C1/3) « SCXSB 00000180 

T1C2.1) = -SCXCA-CBXSAXCC 00000190 

Tl(2,2) * -SCXSA+CBXCAXCC 00000200 

T1 ( 2, 3) = CCXSB 00000210 

T1C3»1) = SBXSA 00000220 

Tl(3,2) * -SB*CA 00000230 

Tl(3,3) * CB 00000240 

RETURN 00000250 

END 00000260 


39 



ooooo 


**x* TSO FOREGROUND HARDCOPY xxxx 
DSNAME=XRHSC. ANTENNA . FORT 


(TRAN5P ) 


SUBROUTINE TRANSPCX# TZ) 00000010 

CXKXXXXKXMX»MKXM*«MKKKXKXKKKXKXKKXXKXKXXXX«MMXKKXXXXXXXXXXMXMfc»l>»MKXXXKXOOOOD02O 


00000030 

SUBROUTINE TRANSP . USING ROTATION OPERATOR TO GENERATE AN 00000040 

INVERSE MATRIX FOR ANT I -DIRECTIONAL 00000050 

TRANSFORMATION 00000060 

00000070 

C*XXXXXXXXXXXXXXXXXXM««*ll»»»*>«MKX««XKMKXMM»«ftXft«K» »*•««»•!••• KM MMKXKXXXXX00000 OS 0 


IMPLICIT REALMS LA-H,0-Z) 
DIMENSION X(3.3J.T2(3,3) 
T2(l,l>* X(l,l) 

T2( 1 # 2) * X( 2» 1 ) 

T2( 1 , 3) = X( 3» 1 ) 

72(2,1) = X( 1 , 2) 

T2(2,2) * X(2,2) 

T2(Z,3) * X(3»2) 

T2( 3, 1 ) * X( 1 , 3) 

T2C 3, 2) «= X(2,3) 

T2( 3, 3) * X(3, 3) 

RETURN 

END 


00000090 

00000100 

00000110 

00000120 

00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 

00000210 


non 


**** TSO FOREGROUND HARDCOPY xxxx 
DSNAME = XF.HSC . ANTENNA . FORT 


(CROSS ) 


SUBROUTINE CROSS( Y,Z, T3) 00000010 

CXXXXXXXXXXXXXXXXXXXXXXXXXXKXXKXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX* 00 000020 

00000030 

SUBROUTINE CROSS : TO PERFORM MATRIX MULTIPLICATION 00000090 

00000050 


CXXXXXXXKXXXXXXXXXXXXXXXXXKXXXXXXXX»XXKK*XXXXXXXX**X*XX*X*XXXXXXXXXXXXXX00000060 


IMPLICIT REAL*« <A-H,0 -Z) 

DIMENSION Y(3,3),Z(3,3),T3<3.3) 

T3(l» 1 ) = Y(1,1)X2(1,1)+Y(1.2)XZ<2,1)+Y(1,3)XZC3,1) 
T3(l,2) = Y<1,1)XZ(1,2)+Y(1,2)XZC2,2)+Y(1,3)XZ(3,2) 
T3(l,3) » Y(l,l)xZ(l , 3)+Y( 1 *2)*Z(2» 3)+Y(l » 3)*ZC3» 3) 
T3(2,l) » Y(2,nxZ<l.l)+Y<2,2)*Z(2,l)+Y<2,3)XZ<3,l) 
T3(2,2) * Y(2,1)XZ(1,2)+Y<2,2)XZ<2,2)+Y(2,3)XZ(3,2) 
T3(2, 3) * Y(2,1)*Z(1»3)+Y(2*2)XZ(2>3)+Y(2*3)XZ(3»3) 
T3<3,1) » Y(3»l)xZ(l»l)+Y(3,2)*Z(2*l)+Y(3#3)xZ(3#l) 
T3C3.2) « Y(3,1)XZ(1,2)+Y(3.2)KZ<2,2)+Y(3,3)XZ(3,2) 
T3( 3, 3) * Y(3,1)*Z<1,3)+Y(3,2)*Z(2,3)+Y(3,3)*Z(3,3) 
RETURN 
END 


00000070 

00000080 

00000090 

00000100 

00000110 

00000120 

00000130 

ooooouo 

00000150 

00000160 

00000170 

00000180 

00000190 



oooo 


XKXK TSO FOREGROUND HARDCOPY XXXX 
DSNAME-XRHSC . ANTENNA . FORT 


(ROT 


) 


SUBROUTINE ROT ( D , EX , EY . EZ # RX , RY , RZ 1 00000010 

CXKKMKXXXXKKKKXXXXXXKXKKKKHMKXMMMHMMRKKKHMMMXKXKKKKMMNKNKKKKXKKXXKKKKKKX00000020 

00000050 

SUBROUTINE ROT • TO PERFORM ROTATION FROM ONE COORDINATE 00000040 

SYSTEM TO ANOTHER ONE 00000050 

00000060 

CXXXXKXKKKKXMKXMXXXXXXKXXXXKXXXXXXKKKXXXXXXXXXXXXKKXXXKKKKKKKKKXKKXXMKXXOOO 00070 
IMPLICIT REALMS (A-H,0-Z) 00000080 

DIMENSION D<3,3) 00000090 

RX « D(1,I)nEX+D(2,1)*EY+D(3,1)nEZ 00000100 

RY • D(1,2)MEX+D(2,2)*EY+D(3,2)MEZ 00000110 

RZ « D(1,3)XEX+D(2,3)MEY+D(3,3)XEZ 00000120 

RETURN 00000150 

END 00000140 
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CHORAL PAQE tS 

OF POOR Q U/^-fTY 


**## TSO FOREGROUND HARDCOPY MIX 

DSNAME=XRHSC. ANTENNA. FORT (ATXPW ) 


SUBROUTINE ATXPW 00000010 

XXXXKXKXkKKKKKXaKKKXXXKXXXXXKKXKXKXXXXXXXXXKXXXXXXXXKXXXXXXXK«**XaxKXKX00000020 

00000030 

SUBROUTINE ATXPW « TO ROTATE THE PLANE-WAVE INCIDENT ANGLES AND 00000040 

POLARIZATION ANGLE FROM INTIAl TO FENCE 0C000050 

SYSTEM 00000060 


00000C7C 

XXXKXXKXXXKXKXKXKKX XX KM X * XXX X HXXKXXXKKXKXXXK XXXXX VKXXKX X XX KKXXXXKXXXXXX 0 000 0 0£ 0 
IMPLICIT REAL *8 (A-H.O-Z) 00000090 

REALXS KIX, KIY, KIZ , KFX, KFY , KFZ, KHX, KHY , KHZ 00000100 

COMMON /INPUT 2/ ALPHAI , BETAI , DELTAI , AL FI , BET1 , GAM1 00000110 

PIE. SPM, RAD. XK.XLM, PA, PB, PDB 00000120 

RIA(3,3),RIF(3,3),RFA(3,3),RAI(3.3).RFI< l. I 00000130 

KIX, KIY , KIZ, KFX, KFY, KFZ, KHX, KHY, KHZ 000001 A0 

HIX.HIY.HIZ 00000150 

ALPHA, BETA, SINALP,COSALP,SINBET,COSBET C0000160 

SINDEL , C0SDEL 00000170 

DSIN(ALPHAI) 00000180 

DCOS(ALPHAI) 00000190 

DSIN(BETAI) 00000200 

DCOS(BETAI) 00000210 

DSIN(DELTAI) 00000220 

DCOS(DELTAl) 00000230 

CAIXCBI 00000240 

SAIXCBI 00000250 

SBI 00000260 

CALL ROT(RJF,KIX,KIY,KIZ,KFX,KFY,KFZ) 00000270 

BETA = DASIN(KFZ) 000002B0 

ALPHA « DATAN2C KFY, KFX) j 00000290 

SINALP * DSIN(ALPHA) 00000300 

COSALP * DCOS( ALPHA) 0000D510 

SINBET « DSIN(BETA) 00000320 

COSBET = DCOS(BETA) 00000330 

HIX = CDIX(-SAI)+SDIXC-CAIXSBI) 00000340 

HIY * CDIXCAI+SDIXC-SAIXSBI) 00000350 

HIZ = SDIXCBI • "• 00000360 

CALL ROT(RIF,HIX,HIY,HIZ,HFX,HFY,HFZ) 00000370 

SINDEL * HFZ/COSBET 00000380 

COSDEL « (HFY+SINDELXSINALPXSINBET)/COSALP 00000390 

WRITE (6,20) COSDEL, SINDEL 00000400 


COMMON /ANTO/ 
COMMON /ANTROP/ 
COMMON /ANTXP1/ 
COMMON /ANTXP2/ 
COMMON /ANTXP3/ 
COMMON /ANTXP4/ 
SAI 
CAI 
5BI 
CBI 
SDI 
CDI 
KIX 
KIY 
KIZ 


20 FORMAT C • 0 ', 15X, • COSDEL = VFS.S,' SINDEL = *,F8.3) 
RETURN 
END 


00000410 

00000420 

00000A30 
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xxxx TSO FOREGROUND HARDCOPY xxxx 
DSNAME s XRHSC . ANTENNA . FORT 


(ATAPER ) 


SUBROUTINE ATAPER 00000010 

CKXXXXXXKKXKX«KXMMKKXKXKXXXXXXXXXKXKXXXXXXKKKXXXKKXXXXXXKKXXXXKXK««ftK»KK00000020 

00000030 

TO SUBDIVIDE THE REFLECTOR APERTURE INTO 
SMALL DIFFERENTIAL AREAS AND EVALUATE 
THE COORDINATE, UNIT NORMAL, AND DIFFEREN- 
TIAL AREA FOR EACH SUBDIVISION AT APERTURE 
SURFACE 


SUBROUTINE ATAPER 


CXXXXXXXXXXXXXKMKXKXXXXKKXXXKXXKXKKXKjMtXXXKXXXXKKXXKMKMXXXXXXttftHXI 

COMMON /INPUT1/ SIGMAM,SIGMA075nX FR, F 
COMMON /ANTO/ PIE, SPM, RAD,XK,XLM, PA, PB, PDB 
COMMON /ANTPA1/ XL C 11993) ,YIC11993) , ZL C 1 1993) , INDEX 
COMMON /ANTPA2/ XNIC 11993) ,XNJC11993) , XNKC 11993) , DSSC 11993) 
DIMENSION SZEPRC 1 1993 ), CZEPRC 11993) 

DSIGO * XLMxXLI 
SIGMAP > SIGMAO 

XMI1 « CSIGMAM-SIGMAO)/DSIGO+l . 0 
IXMII * XMII 

DSIG » CSIGMAM-SIGMAOl/IXMII 
INDEX * 1 

DO 500 IPB «1, IXMII 

SIGMA « SIGMAP+DSIG 

SIGMR > SIGMA-0. 5XDSIG 

XMI « 2.0XPIEXSIGMA/DSIG/A.O+1.0 

JXM * XMI 

IXMI * JXMXA 

JJXM * 2KJXM 

JJJM « 4XJXM 

DZETA « PIE/ JJXM 

ZEPR «-PIE/JJJM 

IN * 0 

DO 200 INN=1 , JXM 
IN « IN+1 • 

ZEPR « ZEPR+DZETA 
SZEPRC IN) *DSINCZEPR) 

CZEPRC IN) *DCOS(ZEPR) 

200 CONTINUE 
JX2 « 2XJXM 
JXM2 * JXM+1 

NK. = 0 

DO 210 IN2 S JXM2, JX2 


SZEPRC IN-NK) 
-CZEPRC IN-NK) 


SZEPRC I N2 ) 

C2 EPRi I N2 ) 

NK « NK+1 
210 CONTINUE 
IN * 0 
JX3 « 5XJXM 
JXM3 * 2XJXM41 
DO 220 1N3 * JXM3, JX3 
IN « IN+1 

SZEPRCIN3) • -SZEPRC IN) 
CZEPRC INS) * -CZEPRC IN) 
220 CONTINUE 

JX4 «= 4XJXM 
JXM4 « 3XJXM+1 


00000040 

00000050 

00000060 

00000070 

00000080 

00000090 

(X00000100 

00000110 

00000120 

00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 

00000210 

00000220 

00000230 

00000240 

00000250 

00000260 

00000270 

00000280 

00000290 

00000300 

00000310 

00000320 

00000330 

00000340 

00000350 

00000360 

00000370 

00000380 

00000390 

00000400 

00000410 

00000420 

00000430 

00000440 

00000450 

00000460 

00000470 

00000480 

00000490 

00000500 

00000510 

00000520 

00000530 

00000540 

00000550 

00000560 


44 


NJ * 0 

DO 230 IN9*JXM9» JX9 
SZEPR(IN9) « -SZEPRCIH-NJ) 

CZEPRHN9) * CZEPR(IN-NJ) 

NJ * NJ+1 
230 CONTINUE 

DO 300 JPB * I , JJJM 

XL( INDEX) « SIGMR*SZEPR( JPB) 
YUINDEX) « -SIGMRXCZEPR(JPB) 

ZL( INDEX) « 0.0 

PAR ■ DSQRT(SIGMRmk2+9 . 0XFXX2) 

XNI( INDEX) « -XUINDEX)/PAR 
XNJ( INDEX) « -YL (I NDEX )/PAR 
XNK< INDEX) « 2 . ONF/PAR 
DS « PIE*<SIGMA*«2-SIGMAPXX2)/IXMI 
DSS( INDEX) • DS 
INDEX » INDEX* 1 
300 CONTINUE 

SIGMAP • SIGMA 
500 CONTINUE 

INDEX • INDEX - 1 
WRITE (6,600) INDEX 
600 FORMAT < • 1 ' , 20X, • INDEX * ’,110) 

RETURN 

END 


00000570 

000005S0 

00000590 
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*XXX TSO FOREGROUND HARDCOPY XXXX 
DSNAME*XRHSC . ANTENNA . FORT 


(ATRAT ) 


SUBROUTINE ATRAT 00000010 

XKMNKXNNXXRXKKXXKXKXXKXXXXXXXXXXKXKKXXXKXKKKKKKKXKKXXKXXXKXKKKKKKXO 00 000 ^° 

00000030 

SUBROUTINE ATRAT i TO TRANSFORM THE POINT IN THE CARTESIAN 00000060 

COORDINATE SYSTEM AT REFLECTOR APERTURE 00000050 

TO THE POINT IN THE SPHERICAL COORDINATE 00000060 

SYSTEM AT FENCE 00000070 


CXXXXXXXXXXXRRXXXXXNNXXXXXXKXXKXXXXXXXKXXXXKXXXXXXXXXXXXXKXKXXKXRX 
IMPLICIT REALX8 (A-H.O-Z) 

COMMON / 1 NPtfT3?*~~A L F2rB ETC; G A ME » *T A > ~YT A , Z T A ,TT F , VT F , ZT F 
COMMON /ANTROP/ RIA(3,3),RIF(3,3),RFA(3,3),RAI(3,3),RFIC3,3) 
COMMON /ANTPA1/ XL( 11993), YL( 11993), ZL< 11993), INDEX 
COMMON /ANTRAT/ RDPR,ZDPR, THETA 
COMMON /ANTLP1/ LOOP1 
C WRITE (6,10) LOOP1 
10 FORMAT ( *0 * , IX, 15) 

TAX « XL(LOOPl) 

TAY « YL(LOOPl) 

TAZ « ZLCLOOP1) 

CALL ROTCRAI, TAX, TAY, TAZ, TIX, TIY, TIZ) 

TIX « TIX+XTA-XTF 
TIY « TIY+YTA-YTF 
TIZ » TIZ+ZTA-ZTF 

CALL ROT(RIF,TIX,TIY,TIZ»TFX,TFY,TFZ) 

RDPR ■ DSGRT(TFXXX2+TFYXX2) 

ZDPR ■ TFZ 

THETA * DATAN2(TFY,TFX) 

RETURN J 

END i 
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oooo 


ORiGiMAL PAQ£ tS 

OF POOR QUALITY 


KXKK TSO FOREGROUND HARDCOPY XRXX 

DSNAME'XRHSC. ANTENNA. FORT ( ATPL ) 


SUBROUTINE ATPl 


00000010 


CXKXXXXXXKKKXRKKKXKRKHXKXXKXXXXXXKXXXXXXXXKXXKXXXXXXXXXXXKXKKKXXKKXKKKMX00000020 


00000030 

SUBROUTINE ATPL . TO COMPUTE THE PLANE-WAVE GCP) AND G(Q) 00000040 

VALUES AS PART OF SOMMERFELD SOLUTION 00000050 


CXKXXXKKXKXXXXXXXXXXXXKXKXXKXKKKKKKXRHKKRKKXXRRRRRRRKRRXRKKK 
IMPLICIT REAL KB (A-H.O-2) 

COMMON /ANTO/ PIE,$PM,RAD,XIC,XLM,PA,PB,PDB 
COMMON /ANTXP3/ ALPHA, BETA, SINALP,COSALP,SINBET,COSBET 
COMMON /ANTRAT/ RDPR,»PR, THETA ~ 

COMMON /ANTPL/ GPR, OPI , GQR, GQI 
DIMENSION GR(2),GI(2) 

TWOPI « 2.0XPIE 
ARG1 • PIE/ A . 0 
C1R * DCOS(ARGl) 

C1I * DSIN(ARGl) 

AR1 « CTHETA+ALPHA)/2.0 
AR2 « CTHETA-ALPHA)/2.0 
Q * -DSQRT(2.0XXKXRDPRXC0SBET)XDC0S(AR1) 

P ■ -DSORTC2.0XXKXRDPRXCOSBET)XDCOS(AR2) 

DO 50 I « 1,2 
U ■ 1.0 

IF C-Q.LE.0.0) U « 0.0 
SGN *1.0 

IF (Q.LT.0.0) SON « -1.0 

ARG2 > QXX2 

AG2TP * ARG2/TW0PI 

IAG2TP « AG2TP -4- . 

RMD • AG2TP - IAG2TP * 

AG2MD > RMDXTWOPI -i 

C2R « DC0SCAG2MD) 

C2I « -DSINCAG2MD) 

SOPIE « DSQRTCPIE) 

SQPIE2 * DSt)RT(PIE/2. 0) 

SQ2PIE * DSQRTL2. OXPIE) 

C3R « S0PIEXUX(C2RXC1R-C2IXC1I) 

C3I * SQPIEXUXCC2RXC1I+C2IXCIR) 

CALL ATCS ( C , S , ARG2 , AG2MD ) 

CAR * SQPIE2/2.-S0PIE2XC 
C4I r SQPIE2/2.-S0PIE2XS 
C5R « SGNXCC2RRC4R-C2IXC4I) 

C5I * SGN»(C2RRC4I+C2I*C4R) 

GRC I ) * C3R+C5R 
GICI) = C3I+C5I 

Q a p 

50 CONTINUE 
GOR * GR(1) 

GDI « GI ( 1 ) 

GPR « GRC 2) 

GPI « GIC2) 

RETURN 

END 
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**** TSO FOREGROUND HARDCOPY xxxx 
DSNAME'XKHSC. ANTENNA. FORT 


(ATCS 


) 


SUBROUTINE ATCS(C, S, X, XMD ) 

IMPLICIT REAIX8 <A-H,0-2) 

Z*ABS(X) 

ZMD«ABS(XMD) 

2 IFCZ-4.) 3,3, A 

3 C*DSQRT(Z) 

S»Z*C 

Z«Z*Z 

ZMO>ZMDXZMD 

C*C*( (((((. 50998 348 E-l 0 «Z- . 10140729E-7)xZ+. 11605284E-5)*Z 
I ~ • 85224622E-4 ) «Z-* . 3fc938586E-2)XZ-.079788405)XZ+. 79788455) 
S=S*(<((C-. 6677744? E-9«Z*.112Z5331E-6)XZ-.10525853E-4)KZ 
1+. 6 043537 1E-3)XZ- If. 997 11 OF. -1 )*Z+ . 26596149) 

RETURN 

A D*DCOS(ZMD> 

S*DSIN(ZMD) 

Z«4./Z 


A*((((C(< .876 8258 3E-3XZ- . 41 692894 E-2)*Z* . 797 09430 E-2)*Z- 
i •«I 9 2 8<, 1 1 E-2J*Z».30953A12E-3)XZ+.59721508E-2)KZ-.1606A281E-A)KZ- 
2 . 024933215>*Z- . 44440909E-8 


B=< (<( (< -.66 339256 E-3KZ+. 34 014090E-2)xZ- . 72716901E-2)*Z+ 

1 .7 4282459E-2)XZ-.40271450E-3)XZ-. 931491 05E-2)*Z-.12079984E-5)*Z+ 
2 . 1994711 
Z=DSQRT(Z) 

C'.5+Z«<DKA+S*B) 

S* . 5+ZXCSKA-DxB) 

RETURN 

END 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


xxxx TSO FOREGROUND HARDCOPY XXXX 

DSNAME*XP.HSC. ANTENNA. FORT (ATCY ) 


SUBROUTINE ATCY 
CXXXXXXXKXXXKMXKKXXXXKX: 


00000010 
XXXKKXXXXXKO 0000 020 


SUBROUTINE ATCY i TO EVALUATE H-PLANE SOMMERFELD SOtUTION 

Ckxkxkkxxxxxxkkxxkxkkxkmkmmxxkkkkkkxkxxxxxxxxxxkkkkxxxxx* a ■ *«»« 
IMPLICIT REALXB (A-H.O-Z) 

COMMON /INPUTX/ ANGA, ANGB, DB, NANG, XDF 

COMMON /ANTO/ PIE,SPM,RAD,XK,XLM,PA,PB,PDB 

COMMON /ANTXPJ/ ALPHA, BETA, SINALP,COSALP, SI NBET , COSBET 

COMMON /ANTXP4/ SINDEL,COSDEL 

COMMON /ANTRAT/ RDPR,ZDPR, THETA 

COMMON /ANTPL/ GPR, GPI , GOR, GOI 

COMMON /ANTCY/ FHXR, FHXI, FHYR,FHYI , FHZR, FHZI 

ARC • PIE/A. 0 

C1R • DCOS(-ARG) 

Cll » DSIN(-ARG) 

ARG3 • XKX(RDPRXCOSBET-ZDPRXSINBET) 

C5R • DCOS(ARG3> 

C5I « DSINUR03) 

C6R « CC5R*ClR-C5lxClIVDS9RTCPIE) 

C6I • (C5R*C1I+C5IXC1R)/DSQRT(PIE) 

IF (XDF.E9.1.0) THEN 

C7TMP « DSQRT ( 2 . 0/(XK*RDPRxcOSBET ) ) 

ELSE 

C7TMP » 0.0 
END IF 

SAH « DSIN( ALPHA/2.0) 

CAH « DCOSCALPHA/2. 0) 

STH « DSINLTHET A/2.0) 

CTH • DCOSCTHET A/2.0) 

C7ISC * C7TMPNSAHXCTH 
C7ISS * C7TMPXSAHXSTH 
C7ICS « C7TMPXCAHXSTH 
C7ICC « C7TMPXCAHXCTH 
C8R * SINALPXCGPR+GQR) 

C8I « SINALPX(GPI*GQI)+C7ISC 
HXER » -( C6RXC8R-C6I»C8I ) 

HXEI « -( C6RKC8I+C6 IXC8R ) 

HXE « HXERXX2+HXEI**2 
C9R * COSALPX(GPR-GQR) 

C9I » C0SALPXCGPI-G0D-C7ISS 
HYER « C6RXC9R-C6IXC9I 
HYEI * C6R«C9I+C6I*C9R 
HYE * HY£RX»2-»HYElx«2 
HZER « 0.0 
HZEI -0.0 

HE « DSQRT C HXE+HYE+HZE) 

Cl OR • CO$ALP*(GPR+GQR)XSINBET 

C10I * (COSALP*(GPI+GQI)+C7ICC>*SINBET 

HXHR * -< C6RXC10R-C6IXC10I ) 

HXHI * -( C6RXC10I+C6IXC10R ) 

HXH * DSQRT ( HXHRX*2+HXHI**2 ) 

C11R « SINALPX(GPR-GQR)XSINBET 

C11I * (SINALPX(GPI-GQI)+C7ICS)XSINBET 

HYHR • -( C6RXC11R-C6IXC11I ) 
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HYHI * 
HYH * 
C12R • 
C12I ' 
HZHR * 
HZHI * 
HZH « 
FHXR * 
FHXI * 
FHYR * 
FHYI « 
FHZR * 
FHZI * 
RETURN 
END 


-< C6RkC11I+C6I*CUR ) 

00000570 

DSQRT(HYHRK«2*HYHI«*2) 

00000580 

(Gpr+oqrjucosbet 

00000590 

(GPI+CQI )*C0SBET 

00000800 

C.6R*C12R-C6I*C12I 

00000610 

C6R*C12I->C6I*C12R 

00000620 

DSQRT(HZHRX*24HZHIX*2) 

00000650 

SINDELXHXHR4C0SDELKKXER 

00000640 

SINDEL*HXHI+COSD£L*HXEI 

00000650 

SINDEIKHYHR+COSDELXHYER 

00000660 

SINDELKHYHI+COSDELXHYEI 

00000670 

SINDELXKZHR4COSDELKHZER 

00000680 

SINDEl*HZHI+COSDEl*HZEI 

00000690 


00000700 


00000710 


50 


noon 


xxxx TSO FOREGROUND HARDCOPY XXXX 
DSNAME*XRHSC .ANTENNA . FORT 


(ATRFA ) 


SUBROUTINE ATRFA 00000010 

CXKKKKXKXKKXKXXXXXXXKXXXKKXKXKXKKKKKXXKXKXXKXXXKK«MXKKKKKKKXXXXXXXKXXKXK0000C020 

00000050 

SUBROUTINE ATRFA . TO ROTATE H-PLANE SOMMERFELD SOLUTION 

FROM FENCE TO REFLECTOR APERTURE 


CKXXKKKKXXXXKKKXXXXXKKKMKXKKKKKKKXKKKKKKKKXKKKXXXXKKXXKMKXRMXMMKXX1 

IMPLICIT REALXB (A-H,0-Z) 

COMMON /ANTROP/ RIAC3,3),RIF(3.3),RFA<3,3),RAIC3,3),RFI(3.3> 
COMMON XANTCY/ FHXR» FHXI , FHYR, FHYI , FHZR » FHZI 

COMMON /ANTRFA/ AHXR , AHXI > AHYR , AHYI , AHZR , AHZI 

COMMON XANTRF1/ HXRC11993),HYR(11993),HZRC11993) 

COMMON /ANTRF2/ HXI (11993) , HYK11993) ,HZI( 11993) 

COMMON ✓ANTLPl/ LOOP1 

CALL ROT ( RFA , FHXR , FHYR . FHZR , AHXR , AHYR , AHZR ) 

CALL ROTC RFA, FHXI . FHYI , FHZI , AHXI , AHYI , AHZI ) 


HXR(LOOPl) 

HYR(LOOPl) 

HZR(LOOPl) 

HXKL00P1) 

HYKLOOP1) 

HZI(LOOPl) 

RETURN 

END 
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ORIGINAL PAGE IS 

OF POOR QUALITY 


***« TSO FOREGROUND HARDCOPY XXXX 

DSNAME=XRHSC. ANTENNA. FORT (ATCORR ) 


SUBROUTINE ATCORR 00000010 

Cxxxxxxxxxxxxxxx x*xx*» »¥»»»» ••*«x#*xxx««x*xxxxx**xxxxxxxxxxxxxxxxxxxxx **00000020 

00000030 

SUBROUTINE Arc OR* ' f- cut ATF i ME CCRPf l AT] CN OF THE 00000040 

r.jiif wvt function and the incoming 00000050 

ri F *RaC.T EE HfCTROMAC.Nt.TIC FIELD AT 00000060 

ANIENNA APERTURE, THE TRANSMJTTFD 00000070 

REGION. OOOOOOSO 

00000090 


CXXXXXXXKKXXXXKKKXXXXXXXMKXXXXKKXXXXXXXXXXXXKXKKXXXXXXKXXXX*R«y««"X**»«M00000100 


IMPLICIT REALX8 (A-H,0-Z) 

COMMON /INPUTX/ ANGA, ANGB, DB,NANG,XDF 
COMMON /ANTO/ P I E, SPM, RAD, XK, XLM, PA, PB, PDB 
COMMON /ANTPA1/ XL C 1 1993 ), U ( 1 1993) , ZL C 1 1993 ), INDEX 
COMMON /ANTPA2/ XNI ( 1 1 993) , XNJC 1 1993) , XNKC 1 1993) , DSSC 1 1993 ) 
COMMON /ANTRF1/ HXR( 1 1993) , HYR( 1 1993) , HZR( 1 1 993) 

COMMON /ANTRF2/ HXI ( 11993) , HYI( 11993) , HZI( 11993) 

COMMON /ANTPW1/ HXARC 1 1 993) , HXAI C 1 1993) , HYAR< 1 1993) 

COMMON /ANTPW2/ HYAI(11993),HZAR(11993),HZAI(11993) 

VOCR * 0.0 

VOCI * 0.0 

DO 800 I = 1, INDEX 

RDNDSS « DSSC I) 

FNR * HXR(I)XHXARCI) - HXI(I)XHXAICI) + HYRCI)XHYAR(I) 

1 -HYI (I)XHYAI(I) + HZRCI)XHZAR(I) - HZI ( I )XHZAI ( I ) 

FNI * HXR(I)XHXAICI) ♦ HXI(I)xHXARCI) + HYR(I)xHYAICI) 

2 +HYI(I)*HYAR( I) +HZRCI, ZAICI) + HZI(I)XHZARCI) 

DVOCR > FNR X RDNDSS 

DVOCI * FNI X RDNDSS j 

VOCR * VOCR + DVOCR 
VOCI * VOCI ♦ DVOCI 
800 CONTINUE 

VOC «= DSORTC VOCRXVOCR + VOCIXVOCI) 

VOCL = 20.0XLOG10CVOC) 

WRITE (6,850) J> B/RAD, VOC, VOCL 

350 FORMAT ( , 0 , ,20X,*PB = '^8.3,' VOC = ',E12.<, ( VOCL = ’,E12.6) 
RETURN 
END 
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xxxx TSO FOREGROUND HARDCOPY xxxx 
DSNAME s XRHSC. ANTENNA . FORT 


(ATPWH ) 


SUBROUTINE ATPWH 00000010 

XXXXXXXXXX XXXX XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XX xxxxxxo 0000 0?0 

00000030 

SUBROUTINE ATPWH . TO EVALUATE THE SHEET-CURRENT AT EVERY 00000040 

POINT ON THE REFLECTOR APERTURE 00000050 

00000060 


CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX00000070 


IMPLICIT REALMS (A-H.O-Z) 

COMMON /ANTO/ PIE,SPM,RAD,XK,XLM,PA,PB,PDB 
COMMON /ANTPA1/ XL ( 11993) ,YL( 11993) ,ZL( 11993), INDEX 
COMMON /ANTPA2/ XNI ( 11993) , XNJC11993) ,XNK(11993) , DSSC 11993) 
COMMON /ANTPW1/ HXAR( 11993) , HXAIC 11993) , HYAR( 11993) 

COMMON /ANTPW2/ HYAI(11993).HZAR(11993),HZAiai993) 

CPA = DCOS(PA) 

SPA « DSIN(PA) 

CPB = DCOSCPB) 

SPB « DSINCPB) 

CACB * CPA X CPB 
CASB x CPA X SPB 
SACB x SPA X CPB 
SASB * SPA X SPB 
DO 300 LP* 1, INDEX 

PS « XL(LP)*CACB + YL(LP)*SACB + ZL(LP)xSPB 
PKS x XK * PS 
CKS • DCOS(PKS) 

SKS « DSIN(PKS) 

HXAR(LP) « -SPA x CKS 
HYAR(LP) « CPA * CKS 
HZAR(LP) « 0.0 X CKS 
HXAI(LP) « SPA X SKS 
HYAI(LP) « -CPA X SKS 
HZAI(LP) * 0.0 X SKS 
300 CONTINUE 

RETURN * 

END 


00000080 

00000090 

00000100 

00000110 

00000120 

00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 

00000210 

00000220 

00000230 

00000240 

00000250 

00000260 

00000270 

00000280 

00000290 

00000300 

00000310 

00000320 

00000330 

00004340 

00000350 


i 


53 



